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Abstract 



Experimental results for condensation in compact heat exchangers show that the heat 
transfer due to condensation is significantly better compared to classical heat exchangers, 
especially when using R134a instead of water as the refrigerant. This suggests that surface 
tension plays a role. Using generalized dimensional analysis we derive reduced model equa- 
tions and jump conditions for condensation in a vertical tube with cylindrical cross section. 
Based on this model we derive a single ordinary differential equation for the thickness of 
the condensate film as function of the tube axis. Our model agrees well with commonly 
used models from existing literature. It is based on the physical dimensions of the problem 
and has greater geometrical flexibility. 

1 Introduction 

Condensation is important in the refrigeration, automotive and process industries, [Caj . Higher 
energy efficiency requirements and the move to more environmentally friendly refrigerants 
increased the need for highly efficient heat transfer for in-tube condensation (and evaporation) 
processes, [Sa]. Improved heat transfer technologies are currently used to save not only energy 
but rather to save space. Over the last decades experimental studies show that the heat transfer 
is better in compact heat exchangers than in classical tube condensers, which made compact 
heat exchangers popular, |Coj . 

The hydrodynamic flow channels in such condensers have diameters in the millimeter range 
and are often inclined to the vertical. The fundamental mechanisms of heat and mass transfer 
as well as of two phase flow in these small channels are not well understood. In this paper we 
study mainly condensation in a vertical tube and make some comments on inclined tubes. 

Literature on condensation in tubes (or channels) with small diameters is on vertical, inclined, 
or horizontal tubes (or channels) and the effect of surface tension is either taken into account 
or not. 

Condensation in vertical tubes is investigated by the following authors, where only the first 
two authors considered surface tension. According to |WaDu2) small surface waves enhance 
the heat transfer mainly due to film thinning effect. [Zh] investigated condensation in vertical 
triangular channels with a diameter between 0.2 mm and 0.3 mm. |Paj showed that the effect 
of interfacial shear stress on the heat transfer depend on the vapor velocity and on the mass 
transfer. [Panj showed that turbulent flow enhances the heat transfer (in a tube with a diameter 
of 24 mm). 
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Some studies on inclined tubes are related to our study. [FiRoj derived an analytical solution 
for condensation in and on elliptical cylinders. He neglected variations of the film thickness 
with the tube radius. [Mo] studied interfacial shear stress without considering surface tension. 
|WaDu] compared heat transfer in horizontal and inclined tubes with constant diameter by 
energy considerations. 

[Sij studied condensation in horizontal parallel plate channels without considering surface ten- 
sion. [WaHo| investigated horizontal micro-fin tubes by dividing the flow in two flow regimes. 
Heat exchangers often operate at moderate Reynolds numbers {Re < 100), with a film thick- 
ness in the millimeter range and small surface waves with wave length in the centimeter range. 
It is well known from linear stability analysis of thin films, that surface tension has a stabi- 
lizing effect on the film, [Be], jYI], and that condensation (opposed to evaporation) tends to 
stabilize the film, |Baj , [MaLe] , [UnTh] , |Sp| . Nonlinear stability analysis for flows along plates 
has shown that stability of the film depends on the frequency of the initial perturbation, see 
[BuBaPa] . [JoDaBa] (evaporation), and [HwWe] (condensation). 

In most of these studies simplified model equations were used. Interestingly, sometimes differ- 
ent equations for the interface between the two phases (jump conditions) are used. Often the 
model equations are simplified intuitively or the method of simplification remains unclear. In 
this paper we derive effective model equations for condensation in vertical tubes with cylindri- 
cal cross sections where surface tension is taken into account. We simplify the model equations 
using generalized dimensional analysis (GDA), an extension of dimensional analysis that allows 
to assign different length scales for each spatial variables. This method has mathematical rigor 
and is very algorithmic, but not well documented in literature, so we review it in this study. 

The rest of the paper is organized as follows: In section [2] the equations for the condensate and 
vapor phase and for the interface between the two phases are given. In section [3] we review 
generalized dimensional analysis (GDA). In section |4] and [5] we obtain reduced model equations 
for the condensate flow, the vapor phase and the jump conditions between the two phases, 
based on the ratio of the film thickness and the tube length, e = H/L, and we obtain the 
nondimensional numbers of the problem. In section [6] we evaluate the dimensionless numbers 
for water and R134a and obtain a single ordinary differential equation for the film thickness 
and in section [7] we compute the heat transfer and compare our results to commonly used 
models from the literature. 

2 Mathematical model 

We consider a vertical tube with cylindrical cross section so that gravity acts in the direction of 
the tube axis and we assume laminar condensate flow. As the velocity of the condensate flow in- 
creases the interface between condensate and vapor becomes wavy with at first two-dimensional 
surface waves, for small surface waves the condensate flow is still laminar, [AlNaPo] . [YoNoNa] . 
Based on this we can assume rotational symmetry so that the variables do not depend on the 
rotation angle. 

In cylindrical coordinates the domain of the condensate flow is given by 

^cond = {(r(t)cosi?,r(t)sini?,z) G | r{t) G [h{z,t),R), -d G [0,27r), z G (0,L)} , 
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Figure 1: Tube with rotational symmetry 



and the domain of the vapor is given by 

^vapor = {(r(t)cosi?,r(t)sini?,z) G | r{t) G {0,h{z,t)), i? G [0,27r), z G (0,L)} , 

where h{z, t) is the thickness of the vapor phase, R is the tube diameter and L is the tube 
length. We denote the thickness of the condensate film by H{z,t) = R — h{z,t). 



Assuming that the condensate and the vapor can be modeled as an incompressible Newtonian 
fluid the mass, momentum, and energy balance equations for the condensate flow in 
cylindrical coordinates are given by 
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where v = [vrir, z, t),Vz (r, z, t)] is the condensate velocity, p{r, z, t) is the condensate pressure, 
r(r, z, t) is the condensate temperature, and g is the gravity. We assume constant density p, 
dynamic viscosity fj,, specific heat capacity c, and thermal conductivity A of the condensate. 



The velocity of the vapor phase is assumed to be small, so that shear stress exerted by the 
vapor on the condensate film can be neglected. The vapor temperature shall be constant, i.e. 
we assume zero heat flux in the vapor phase. The pressure in the vapor phase is determined 
mainly by the hydrostatic pressure, 

-g^ = Pv9, (5) 
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[Fi] . Under this conditions the problem reduce to a single-phase problem and only the equa- 
tions for the condensate phase have to be solved. In the following we refer to equations ([T|) - 
([5]) as the bulk equations. 

The boundary conditions for the momentum equations are: non-homogeneous Dirichlet 
boundary condition at the tube inlet, 



v,{r,z = 0,t)=U{r) , 



(6) 



zero velocity at the tube wall, 

Vr{r = R,z,t) = , (7) 
and the outflow condition at the tube outlet, 



v,ir = R,z,t) = 



(9) 



dvz 
dz 



{r,z = L,t) =0 



(8) 



(10) 



The boundary conditions for the energy equation are: the temperature at the tube wall is 
known, 

T{r = R,z,t)=T^ , (11) 



zero temperature gradient normal to the tube inlet and outlet. 



dT 

dz 



ir,z = 0,t) = 



(12) 



dT 
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{r,z = L,t)=0 



(13) 



and we assume thermodynamic equilibrium, i.e. the temperature at the interface is the satu- 
ration temperature of the vapor. 



T{r = h{z,t),z,t)=Ts . 
Pressure boundary conditions will be discussed in section HI 



(14) 



The mass, momentum, and energy jump conditions for the interface between condensate 
and vapor are with the assumption of no-slip at the interface, ' * =0, constant surface 
tension 7, and after neglecting kinetic and viscous terms in the energy jump condition given as 

[ri - It] • n| = , (15) 
n-Dn =2H-f, (16) 



mv ■ n\\ + 



P 



t-Dn 



q ■ n 



0, 
0, 



(17) 
(18) 



where u is the velocity of the interface, m = p [v — u] ■ n is the volume specific mass flux, 
D = jjL [Vi» -|- (Vt))"^] is the deviatoric stress tensor for an incompressible Newtonian fluid, 
H = — • n is the mean curvature (V5 is the surface gradient), A/i^, is the latent heat of 
vaporization, q = AVT is the heat flux vector, and n, t are the normal and tangent vector 
respectively. For a more detailed discussion of jump conditions see [Dz] . 
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Parameterizing the interface implicitly by F{r, z,t) = r — h{z, t) = we obtain from (jlSp - 

m 




3 Generalized dimensional analysis 

In classical dimensional analysis all spatial variables are scaled by the same length scale. Gen- 
eralized dimensional analysis is an extension of dimensional analysis which allows different 
length scales for each spatial variable, [NeJ. 

Length, mass and time (LTM) are fundamental dimensions of the MKS unit system and the 
c.g.s unit system. Meter, kilogram, second (MKS) are units and centimeter, gram, second 
(c.g.s) are other units. Another set of fundamental dimensions are length, force and time 
(LFT). An equation in which the units balance on both sides of the equal sign is called coher- 
ent. An equation in which the dimensions are equal on both sides of the equal sign is called 
homogeneous^ 

According to Buckingham's Il-Theorem every physically meaningful equation /(oi, . . . , a„) = 
with n variables, where the n variables are expressed in terms of r fundamental dimensions, 
can be rewritten as an equation oi n — r = k dimensionless variables constructed from the 
original variables, -F(7ri,7r2, . . . ,TTk) = 0. Physically meaningful equations have to be invariant 
under a change of system of units and this is used in a dimensional analysis. 

For a generalized dimensional analysis the variables in the differential equations and in the 
boundary conditions are substituted by their product of value and dimension. Because every 
equation must be of dimensional homogeneity this result in equations for the dimensions only. 
The derivatives are not carried out over the dimensions. For example, from the continuity 

^For example, the equation 3 m + 3 cm = 3 m + 3 x 0.01 m = 3.03 m is homogeneous but not coherent. 
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equation, 



rr orr ozz 



we get the normalized dimension equation f ^ 



Z Vr V, 



1. 



A system of physically meaningful ordinary or partial differential equations with n variables, 
where the n variables are expressed in terms of r fundamental dimensions, result in / normalized 
dimension equations 
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with k = ranks. Prom the I coherent equations k are linear independent, so that k nondi- 
mensional variables can be constructed. 

The dimension of a variable aj is a function of power monomials of the dimensions of all 
variables, so that we make the following ansatz: 
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where the coefficients yij are unknown. Substituting the dimensions dj in the dimension equa- 
tions (I23p by their power monomials gives for the ith equation 
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SO that the exponents sum up to zero. Doing this for all / unit equations result in the following 
system of linear equations for the exponents of the dimensions: 

By^=0. 

The Buckingham-n Theorem says that from k = rankS linear independent equations, k 
dimensionless variables can be constructed. The number of fundamental dimensions to be 



6 



prescribed are given by the Buckingham II-Theorem as the number of variables minus the 
number of nondimensional variables, n — k = r. 

We collect the r column vectors corresponding to the fundamental dimensions in the matrix 
R and the remaining k column vectors in the matrix K, so that 

B = [K\R] . 

Note that the r fundamental dimensions have to be chosen such that lankK = lankB. Then 
the k nondimensional variables are given by 

vr,=af^af ...at' , j = h...,k. 

If we assign different units (length scales) to the spatial dimensions the equations are not any 
more coherent and the system By^ = becomes inconsistent. This is the interesting case. Then 
we have with GDA a rigorous method to determine the terms that have to be dropped from 
the model equations in order to recover a consistent system and to determine the characteristic 
dimensionless numbers of the problem. For multiscale problems with several small parameters 
the system may have more than one solution. 



4 Reduced equations and dimensionless numbers of condensate 
flow and vapor 

To obtain the normalized dimension equations of the bulk equations we start by writing 
all independent and dependent variables, material constants, and constants that appear in the 
boundary conditions as product of value (hatted variables) and dimension (tilded variables) , 

r = fr , Vr = Vr Vr , P = P P ■, R = RR , U = U U , 

z = zz , Vz= VzVz , M = A M , H = H H , Ty, = fwfy, , 

t = it, P = PP, 9 = 99, L = LL, Tyj = fy,fyj, 
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Substituting the variables in the continuity equation ([T]) by the product of value and dimension 
gives the normalized dimension equation 

r-'^zvrv'^ = 1. (25) 
From the momentum equations ([2]) and ([3]) we obtain 

fi-^v~^ = l, (31) 

f^'^ 2vrV-^ = 1 , (32) 

vIp'^P=1, (33) 

r'^ z-^ y-^ p fi'^ = 1 , (34) 

f = 1 , (35) 

z~'^ Vz p~^ pg~^ = I ■ (36) 
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Prom the energy equation @ we obtain 

r v-^ = 1 , (37) f2 z-^ v^pc~X-^ = 1 , (39) 

f-^ zvrv-^ = 1 , (38) = 1 . (40) 

From the momentum equation of the vapor phase ([5]) we obtain 

(ISSD, (I27j), (I32]),(I38]) are identical, ([26]), ([HII), dST]) are identical, and ([HOj), ([35]), ([lO]) are iden- 
tical. From the four pressure related dimension equations ([28l) . (p9|) . ([33|) . two are linear 
independent, but before we eliminate linear dependent equations we analyze the boundary 
conditions and study its consequences. 

Now we obtain the dimension equations of the boundary conditions. From the non- 
homogeneous Dirichlet boundary condition in terms of values and dimensions, ([6]), 

Vz Vz {rf, z z = 0, £t) = U U{f f) , 

we obtain that the streamwise velocity has the same dimension as the velocity at the inlet, 

Vz = U ■ (41) 

The dimension equation 5 = is homogeneous and provides no information. 

From the remaining velocity boundary conditions at the wall ([7]) , ([8]) and the outflow condition 

(fTO]l we obtain 

f = R, (42) z = L. (43) 

From the energy boundary conditions at the wall (llip and at the interface (I14p together with 
with h{z,t) = R — H{z,t) we obtain, 

f = f^, (44) h = R, (45) R = H, (46) f = f, . (47) 

The inlet and outlet boundary conditions of the energy equation provide no new information. 

Now, we assign tube radius and tube length to the spatial dimensions, R = H = H 

and L = L. However, with this length scales (units) the linear equation system becomes 
inconsistent. To recover a consistent system of dimension equations we compare the order of 
the terms and drop higher order terms. 

i) The second derivatives in streamwise direction are of order = (H/L)'^ smaller than the 
second derivatives in radial direction and should be dropped, so that (i30|l . (f35|l . and (HOj) 
vanish. Note that eliminating (j35p changes ()36p to 

f-2^,p-vr' = 1 (48) 

ii) The pressure related dimension equations (p8|) . (j29|) and ([33]) . ([M]) contradict dimension 
equation ()25p . Setting the radial derivative of the pressure to zero solves the problem. 
Then the dimension equations (j26p - (j29p vanish. 
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iii) Consequently the condensate pressure is a function of streamwise coordinate and time 
only p = p{z,t) and is determined by the hydrostatic pressure of the vapor || = ^ = 
Py g. The radial velocity is one order smaller compared to streamwise velocity. 



Remark 1. We see that i) - iii) are results of our generalized dimensional analysis. However, 
when Prandl derived his famous boundary layer equations he used i) - iii) as initial assump- 
tions, |Sch] . (Prandtl had a great intuition.) 

The boundary conditions also determine the dimensions of vapor thickness and film thickness, 
h = H and D = H. For simplicity we further assign the dimension of inlet velocity and 
saturation temperature, U = U and Tg = Tg. Then the dimension of temperature and wall 
temperature are T^j = Tg and T = T^. 

With the simplifications discussed above and after eliminating linear dependent equations the 
dimension equations form a system of 7 homogeneous linear equations for n = 12 unknown 
dimensions, B x = 0, with k = lankB = 7, so that B has maximum rank. From this equation 
system k = 7 nondimensional variables can be constructed from r = n — k = 5 fundamental 
dimensions, which have to be chosen such that the matrix Rin B = [K\R] fulfills the condition 
ranki^ = rankS. 

Some dimensions are already determined by the boundary conditions. In addition to the spatial 
dimensions and the dimension of the inlet velocity we chose the dimensions of the two material 
properties density and heat capacity as fundamental dimensions, so that 
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Setting p = p, c = c we obtain with r = H, I = L, Vz = U after applying Gaufiian elimina- 
tion to B = [K\R] from the row-reduced echelon form of B the following nondimensional 
variables: 
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we get with /i 
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Using this we obtain from ([1]) ([5]) the reduced mass, momentum, and energy equations 

in dimensionless form as the following system of partial differential equations: 
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and the momentum equation of the vapor flow, 

dpv 
dz 



Pv9 



(49) 
(50) 
(51) 
(52) 

(53) 



These equations are the zero and first order equations for the bulk flow. 

Remark 2. Nusselt's theory of condensation along a flat plate is based on zero order equations, 
|Nu] . Our model is an extension of Nusselt's theory for condensation in a tube of rotational 
symmetry. 



5 Reduced jump conditions and dimensionless numbers of the 
interface equations 

Now we use GDA to derive reduced jump conditions. We have three additional variables, 

h = hh, 7 = 77, Ahv = AhAh. 

First, we obtain the dimension equations of the jump conditions. All jump conditions 
(fT9|) - (122D contain the term \/~^'+~{^)^ ■ With 1 > (^||^ this term reduces to one, so that 

the square root in the jump conditions vanishes. Also the term ^1 — in the tangential 

momentum jump condition vanishes. 

Remark 3. This reflects the assumption of long wavelength approximation. 

Substituting the variables in the jump conditions by the product of value and dimension and 
using the nondimensional numbers vri to 7r4 and ttq we obtain the following jump conditions: 



/ . dh ^ dh\ 
m = pell -Vr + -r^ + fzTTT = peum 
\ dt dz 



(54) 
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2 2^1 ry 

peV m, \ -Vr + Vz-^ j + pU [p - pv) 

2tt2 - \ dVr . f 2dVr , dvA Oh 2 ^Vz ( dh\ . 

-/'^f^ H ^ + + e ^(^^1 I (55) 

7 / od^h 1 

peU-i^ f -e^ - 2 (e^ + |i) + e= = , (56) 

\ or oz \ oz or J oz oz I 

peUAK mAK = -peUcTsX + • (^7) 



The mass jump condition (154p provides no additional dimension equation. 
The first three terms in the normal momentum jump condition (|55p have the dimensions pe^U"^, 
pU'^, and pe^U"^ . Momentum transport due to condensation (or evaporation) and viscous 
normal stress exerted on the interface are of order smaller compared to the pressure terms 
so that from the left hand side only the pressure term remains. From the two curvature terms 
the term related to surface waves is of order two smaller than the term related to the tube 
diameter, so that the first curvature term vanishes. 

The tangential momentum jump condition (j56p reduces to the condition of no shear stress at 
the interface. This yields only a homogeneous dimension equation. 

The second term on the right hand side of the energy jump condition (I57p is of order two 
smaller than the first term and vanishes. 

We obtain two normalized dimension equations, one from the normal momentum jump condi- 
tion and one from the energy jump condition, 

hPvr^ = l, (58) hfr^y^f-^~p\-^KK = l. (59) 

This gives two additional dimensionless numbers from the jump conditions, 

7 
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Defining Weber number and Stefan number as 



We = , St = — , 

7 Ah V 



we get 



Using this we obtain from (j54p - (j57p the reduced mass, normal and tangential momen- 
tum, and energy jump conditions in dimensionless form as the following system of partial 
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differential equations: 



m = peU \-Vr + + V;, — \ , (60) 



dvz 
df 



, (62) 



I I dh ^ dh\ I dT 



The normal and tangential momentum jump conditions are zero order equations. (I62p is the 
condition of no shear stress at the interface. Condensation is a first order process, i.e. mass 
flux terms and phase change terms in the energy jump condition have order epsiloncl 

Remark 4. (j6ip is the Young-Laplace equation. 



6 One ordinary differential equation model 

First, we compare the dimensionless numbers 7r4,7r5,7r6 and vrg, vrg for water and R134a, a re- 
frigerant widely used in automobile air conditioning, to motivate further simplifications of the 
model equations. The dimensionless numbers are the coefficients of the terms in the reduced 
balance equations (f49l) - (1531) and in the reduced jump conditions (i60]l - (f63]l . See appendix lAl 
for material values and constants. The quotient of the two length scales is e = 0.2 10~^. 



dimensionless number 


Water 


R134a 


7r4 = l/(eRe) 


1282 


58 


TTs = l/(eFr) 


8772 


243 


vrg = l/(ePe) 


325 


17 


vrg = 1/We 


1250 


3.37 


vrg = 1/St 


0.002 


0.42 



Clearly gravity forces and viscous forces are dominant compared to inertial and pressure forces 
for both fluids. Conductive heat transport dominates convective heat transport - considerably 
more for water than for R134a. The effect of surface tension is much greater for water than 
for R134a, where the pressure is almost determined by the hydrostatic pressure. At the inter- 
face heat conduction dominates condensation, again the difference is greater for water than for 
R134a. Apparently the flow is the dominating process and the position of the moving surface 
is mainly determined by the solution of the free surface problem. 



^For condensation along a flat plate curvature due to surface waves is the only curvature term and for higher 
Reynolds numbers surface waves may become wavy and cause instabilities. 
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Based on these considerations, we neglect transient and convective terms in the bulk equations 
and the radial velocity in the mass-energy jump condition for the rest of the paper. Then the 
bulk equations reduce to a system of two partial differential equations, 

= ^ - — (64) 
dz ePr ' ^ ' 




= . (65) 
together with p = p{z), i.e. the pressure is the hydrostatic pressure. 

Boundary conditions for velocity and temperature are prescribed at the wall and at the inter- 
face, 

v,{r = R,z) = 0, l-±\{f = h,z) = 0, (66) 

f (f = R,z) = f^, f{r = h,z)=fs. (67) 

Setting = we have from the jump conditions at r = h{z), 

We h 

^ dh _ St df 
^ dz ePe df 

The initial condition for the film thickness h is that the film thickness is prescribed at the inlet, 

/i(r, z = 0) = . (70) 



Differentiating the normal momentum jump condition (168p with respect to z and substituting 
the pressure gradient in ()64p gives 

1 d f ^dvz\ sRe 1 dh Re 

V We" p 5l ~ FV ■ ^ ' 

By integrating (|71|) and (|65|) and evaluating the boundary conditions we obtain the dimension- 
less velocity as 

and the dimensionless temperature gradient as 

df \n{h/R) r' 



where the underlined term in ()72p represents the effect of surface tension and pressure in ^4 
and ()68p . Evaluating (|72p and (I73p at f = h{z) and substituting the results into the mass-energy 
jump condition ([69]) gives 




eRe 1 dh Re \ R^ - h^ h \ dh St f , - f ^ 



2 We dz 2Fr/\ 2 R I dz ePe\n{h/R)h 
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or equivalently 



'dh\ ^ We ^2 5/i _ 2 We St (f, - f^) h 
dzj eFr dz~ e'^ReFe ( + p i^h\\n{h/ R) ' 



Note that with St = ^ and - = 1 - = 4f we have St (f, - f„) = We 

call this new dimensionless number ST. 

This equation is a fully nonlinear first order differential equation for the film thickness. We 
write more conveniently, 

with a = and b = ^a^gp^ • The left hand side can be transformed in an quadratic expression 

(dh h'^V ^.,ru^ 



so that 



= a-±\ a^--bf{h) . 



dz 2 V 4 

Substituting back the parameters a and b we obtain the following quasilinear ordinary differ- 
ential equation: 



dh _We + /We Y 2WeST h 

al ~ elV y ~ Veivy T ~ e^RePe | /^^ In ^ ' 

where the physical relevant solution describing condensation in a vertical tube with small di- 
ameter is the solution with the negative signed root. 

Neglecting the effect of surface tension and pressure in ()74p gives a simpler ordinary differential 
equation, 

dh 2FrST 1 , , 

— = (77) 

dz sRePe (^E^^ + hUniynhh' 

7 Numerical solutions and comparison with other models 

Nusselt's equation for the dimensionless film thickness H oia, condensate flow along a flat plate 
is a function of the dimensionless streamwise coordinate z to the power of one fourth, |Bej| , 
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Water We assume a tube diameter of 7 mm and take water as condensate, with an initial 
film thickness of ho = 0.1 mm. Ah computations are dimensionless. Further details, such as 
wall temperature and the dimensionless numbers of water, can be found in the appendix [XI 
Although the model equations are nonlinear differential equations a standard Runge-Kutta(4) 
method is sufficient. We used Mathematica for the computations. The solutions for the model 




z 



0.2 0.4 0.6 0.8 1 



Figure 2: Water: NuBelt ([78|) . without surface tension ([77|) with surface tension (j76|) 

equations without surface tension and with surface tension are not distinguishable in figure [2j 
The film thickness predicted by (f77l) and (1761) is slightly above the film thickness predicted by 
Nufielt. The reason is that because of the circular tube the condensing mass result in a thicker 
condensate film. 

As a side effect we show in figure [3] that for a diameter of about 60 mm the effect of surface 
tension due to the tube diameter becomes insignificant!! 




z 



0,2 0,4 0,6 0,8 1 

Figure 3: Water: ([78]) , dZZD , ([76]) for d = 60 mm and d = 8 mm 

^Note that the characteristic velocity and length are based on a diameter of 7 mm. 
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R134a For R134a we assume the same tube geometry, initial film thickness, and temperature 
interval as before. The numerical solutions of (j77p and (j76p for R134a are shown in figure HI 
The R134a film is thinner than the water film. Beside this the results are similar as before. The 
difference between Nufielt's solution and the numerical solutions of (j77p and (j76p is smaller, 
because of the thinner film thickness of the R134a film. 




1.5 



0.2 0.4 0.6 0.8 1 

Figure 4: R134a: Nuf^elt ([TSj) . without surface tension ([TTj) . with surface tension ([76|) 



Surface tension We could not observe an effect of surface tension on the film thickness for 
a vertical position of the tube. Surface tension force acts normal to the interface. The effect of 
surface tension can be described by the tendency to minimize curvature (potential energy), and 
the minimal surface of the film is cylindrical. Surface tension is the result of molecular forces. 
In the condensate flow the molecules are surrounded by other molecules and are attracted 
equally in all directions. 

In a vertical tube the film cross sections are always circular, except at the tube inlet where 
the film thickness varies substantial over the tube length. The effect of surface tension to 
minimize curvature results in a more evenly distributed film thickness along the tube length. 
An advantage of our model is that the tube radius does not need to be constant and an analysis 
of the model equations for condensation in a rotational tube where the radius is a function of 
tube length is straight forward. 

Nufielt number, dimensionless heat transfer The heat transfer in the condensate film 
is mostly conductive so that the temperature profile in the film is almost linear and we can 
write 

= aAT . 

Equating g" with the energy jump condition at the interface (where the heat transfer process 
occurs) gives 

a AT = WT ■ n = m A/i„ , 
and in terms of dimensionless variables, 

_ A AT (9f _ Ah^ pUH ^dh 
H dr L dz 
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Multiplying the last equation with gives the local Nufielt number, which is defined as the 
dimensionless temperature gradient at the interface 



_aH _df _ Ah^pUH^ . dh 
~~'dF~ XATL ""dB 



Figure [5] shows the local NuBelt number of Water and R134a as a function of tube length. Note 
that the local NuBelt number is almost inversely proportional to the film thickness. 




Figure 5: NuBelt number of Water and R134a 

Finally we compare our results with the mean NuBelt number Nunu = (3Re)^/'^ for condensa- 
tion along a flat plate according to Nusselt's theory and with the mean Nufielt number given 
by Chen for co-current condensation in tubes [.Bej| , 

NT^chen = (Re-0-4V5.8210-6Re0-8prV3)'^' 

Chen reviewed available experimental information for co-current condensation inside vertical 
tubes. His model is more accurate for Reynolds numbers greater than Re = 30, since it takes 





Water 


R134a 


Nunu 


0.44 


0.15 


Nu 


0.52 


0.23 


Nuchen 


0.74 


0.37 



Table 1: Comparison of mean Nufielt numbers 

into account that at higher Reynolds numbers the film is thinner due to the co-current vapor 
flow. This explains the difference between our model and Chen's model. Our model is based 
on the physics of the problem and has a more general range of applicability. 
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Some remarks on inclined tubes We argued in section [6] that the main forces in the 
condensation process are gravity force and surface tension force and that the position of the 
interface is mainly determined by the solution of the flow problem. 

In an inchned tube gravity force can be divided into two components. One component acts 
in the plane of rotation and the other acts in the direction of the tube axis. Gravity causes 
the condensate to flow in a gathered stream at the bottom of the tube. Surface tension acts 
normally to the surface. Most likely for both fluids the balance of the two forces acting in the 
plane of rotation result in a cross section without rotational symmetry. 

At the surface of the condensate the attraction forces result in a force that is directed inwards 
the condensate film such that the condensate film is pulled into circular shape. Surface ten- 
sion force is the dominant force for water, but is comparably small for R134a, so that for a 
given inclination angle of the tube the cross section of the water condensate film will be more 
cylindrical and the R134a condensate film will mainly fiow at the tube bottom. 
The heat transfer through the condensate is mainly conductive. According to Fourier's law 
the conductive heat flux through a fluid for a given fllm thickness and a constant temperature 
diff'erence is better if the film is thinner. So a better heat transfer can be expected for R134a, 
which would explain experimental results, [Fi|. The model equations for condensation in an 
inclined tube will be partial differential equations. 
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A Material properties and dimensionless numbers 



Water: 

Measured quantities (vapor has saturation temperature): 



vapor temperature 


Ts -- 


= 318.98 




K 


temperature difference 




s 5 




K 


pressure 


p = 


= 10 


103 


Nm-2 


mass flux 


Ml -- 


= 0.05 


10-3 


kgs--'^ 




Mv -- 


= 0.67 


10-3 


kgs-^ 


initial fllm thickness 


Ho p 


a 0.1 


10-3 


m 



Inner tube diameter, characteristic length, characteristic velocity, epsilon: 



di = 2{R - Ho) 

TTdi 



U 



Ml 



E = ^ 



= 7-0.210-3m 

= 0.1 10-3 m 

= 0.023 ms-^ 

_ MMZZiH — n 9 1 n-3 

— 0.5 m — 



Material properties of water at Tg and p: 



density 


Pi 




989.9 






kgm-3 




Pv 




0.068 






kgm-3 


dynamical viscosity 


Pi 




0.60 


10- 


3 


Nsm-^ 




Pv 




162.9 


10- 


3 


Nsm-2 


kinematical viscosity 






0.606 


10- 


6 


9 —1 

m s 








0.011 


10- 


6 


2 —1 

m s 


thermal conductivity 


\l 




0.637 






WK-im-i 








19.98 


10- 


-3 


WK-^m-i 


specific heat capacity 


Cl 




4.179 


103 




Jkg-iR-i 


latent heat of evaporation 






2.393 


103 




Jkg-i 


surface tension 


a 




68.78 


10- 


-3 


Nm-i 



Dimensionless numbers 
Ml 



Re 

Pr 

Pr 
Pe 
St 



LLTrd 
iH 



_ pUH 
P 



~ir 

RePr 
OT _ cp AT 

^-^ - ~AE;r 

We=^ 
Nu = ^ 



3.90 

0.57 

3.94 
15.4 
557 

8.37 
0.7910-3 



R134a (1,1,1,2-Tetrafluorethan): 

Measured quantities (vapor has saturation temperature) 

vapor temperature Ts = 297.15 
wall temperature = 294.15 



K 
K 
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pressure p = 0.65 10^ kgm ^ 

condensate mass flux M; = 0.367 10-^ kgs"^ 

film thickness (vertical tube) Hq = 0.1 .. 0.2 10"^ m 

Characteristic length and characteristic velocity: 

if=4r = 0.1 10-3 m 

U = ^ = 0.142 ms-1 

Material properties at Tg and p according to Tillner-Roth: 



density 


Pi 




1210 




kgm" 


-3 




Pv 




31.39 




kg m" 


-3 


dynamical viscosity 






198.70 


10-6 


kgsm ^ 




P-v 




12.10 


10-6 


kgsm-^ 


kinematical viscosity 






0.164 


10-6 


2 

111 s 


-1 








0.385 


10-6 


m^ s" 


-1 


thermal conductivity 


\i 




82.98 


10-3 


WK- 


-im-1 




Xy 




14.35 


10-3 


WK- 


-1 j^-i 


specific heat capacity 


Cl 




1.421 


103 


Jkg- 


iK-i 




Cy 




1.025 


103 


Jkg- 


■iK-i 


thermal diffusity 


ai 




48.26 


10-9 


2 

111 S 


-1 




ay 




445.9 


10-9 


m^ s" 


-1 


latent heat of evaporation 


hi 




233.1 


103 


Jkg- 


■1 




hv 




411.8 


103 


Jkg- 


1 




Ahy 




178.72 


103 


Jkg- 


1 


surface tension 


a 




8.21 


10-3 


Nm- 


-1 



Dimensionless numbers: 

Re=J%- = ^ = 86.46 

Fr=^ = 20.55 

Pr = ^ = 3.40 

Pe = RePr = 294.2 

St = f^ = 2.36 

ST = = 0.040 

We = P^!/^ = 0.297 



a 
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